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Abstract 

Background: Shigellosis is a major public health concern worldwide, especially in developing countries. It is an 
acute intestinal infection caused by bacteria of the genus Shigella, with a minimum infective dose as low as 10-100 
bacterial cells. Increasing prevalence of Shigella sonnei as the etiologic agent of shigellosis in Malaysia has been 
reported. As there is limited information on the genetic background of S. sonnei in Malaysia, this study aimed to 
characterize Malaysian 5. sonnei and to evaluate the prospect of using multilocus variable-number tandem-repeat 
(VNTR) analysis (MLVA) for subtyping of local S. sonnei. 

Methods: Forty non-repeat clinical strains of S. sonnei isolated during the years 1997-2000, and 2007-2009 were 
studied. The strains were isolated from stools of patients in different hospitals from different regions in Malaysia. 
These epidemiological^ unrelated strains were characterized using biotyping, antimicrobial susceptibility testing, 
pulsed-field gel electrophoresis (PFGE) and MLVA. 

Results: The two biotypes identified in this study were biotype a (n = 29, 73%) and biotype g (n = 1 1, 27%). All the 
40 strains were sensitive to kanamycin, ceftriaxone and ciprofloxacin. Highest resistance rate was observed for 
streptomycin (67.5%), followed by tetracycline (40%) and trimethoprim-sulfamethoxazole (37.5%). All the S. sonnei 
biotype g strains had a core resistance type of streptomycin - trimethoprim-sulfamethoxazole - tetracycline whereas 
the 29 biotype a strains were subtyped into eight resistotypes. All the strains were equally distinguishable by PFGE 
and MLVA. Overall, PFGE analysis indicated that S. sonnei biotype a strains were genetically more diverse than 
biotype g strains. Cluster analysis by MLVA was better in grouping the strains according to biotypes, was reflective 
of the epidemiological information and was equally discriminative as PFGE. 

Conclusions: The S. sonnei strains circulating in Malaysia throughout the period studied were derived from 
different clones given their heterogeneous nature. MLVA based on seven selected VNTR loci was rapid, reproducible 
and highly discriminative and therefore may complement PFGE for routine subtyping of 5. sonnei. 
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Background 

Shigellosis is an acute intestinal infection caused by 
bacteria of the genus Shigella. The main symptom of 
this infection is bloody diarrhoea and the minimum 
infective dose is as low as 10-100 bacterial cells due 
to relative resistance to stomach acid [1]. Shigellosis is 
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a major public health concern worldwide, especially in 
developing countries [2]. The infection is most fre- 
quent in children, the elderly and the immunocom- 
promised [2,3]. A more recent annual estimate of the 
shigellosis burden was at 90 million incidences and 
108,000 deaths [4]. 

Shigella can be differentiated into four species or ser- 
ogroups, S. dysenteriae, S. flexneri, S. boydii, and S. son- 
nei based on biochemical properties and group-specific 
O antigens in the outer membrane of the cell wall 
There is a shift in Shigella dominance from S. flexneri to 
S. sonnei in developing and developed countries [4-6]. In 
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Malaysia, Lee and Puthucheary [7] and Banga Singh et 
al. [8] reported an increasing dominance of S. sonnei as 
the etiologic agent of shigellosis. Apart from the appar- 
ent dominance shift, several studies have described the 
spread of a well-defined pandemic S. sonnei clone since 
1990s which is characterized by common features such 
as biotype g, a particular PFGE Xbal pulsotype, resist- 
ance pattern of streptomycin - trimethoprim- sulfameth- 
oxazole - tetracycline, and the presence of a 2.2 kbp 
class 2 integron [9,10]. 

Bacterial subtyping is frequently applied for outbreak 
investigations and surveillance of infectious diseases. 
Strain- specific fingerprints generated are used to facili- 
tate the identification of disease transmission routes and 
sources [11]. Biotyping is one of the earlier phenotypic 
subtyping methods applied to S. sonnei. The bacteria can 
be subdivided into five biotypes (a, d, e, f, and g) on the 
basis of biochemical properties. This method however is 
not sufficiently discriminative. 

Several genotyping methods with higher discrimin- 
atory power such as pulsed-field gel electrophoresis 
(PFGE) have been applied to subtype S. sonnei [12-14]. 
Since PFGE is a gel-based method, it requires strict ad- 
herence to standardized protocol for reproducible 
results. Standardized PulseNet PFGE protocol has been 
useful for inter-laboratory comparison of results [15]. 
However, due to the monomorphic nature of S. sonnei, 
PFGE occasionally may not be able to distinguish epide- 
miologically unrelated S. sonnei isolates and is not ap- 
propriate for phylogenetic analysis of strains that have 
evolved over a longer time span [16,17]. 

Multilocus variable-number tandem-repeat (VNTR) 
analysis (MLVA) is a highly discriminative sequence- 
based subtyping tool and has been used to study the 
genetic relatedness among S. sonnei. This method is 
based on the inherent variability of short sequences that 
are organized as tandem repeats at multiple VNTR loci. 
The VNTR loci in S. sonnei were found to have different 
degrees of variability [17]. MLVA with four to eight 
highly variable VNTR loci exhibited a discriminatory 
power parallel to or higher than PFGE [18]. Further- 
more, MLVA based on the combination of VNTR loci 
with different variability has also been successfully used 
for phylogenetic analysis of S. sonnei that have evolved 
over different timescales [17]. 

Although S. sonnei is becoming an important etio- 
logic agent of shigellosis in Malaysia, there is limited 
information on the genetic background of local strains. 
Therefore, the objective of the study was to characterize 
local S. sonnei strains by analyzing their biotypes, anti- 
microbial resistance patterns and genotypes. In 
addition, the prospect of using MLVA for routine sub- 
typing of local S. sonnei in comparison with PFGE was 
also evaluated. 



Methods 

Bacterial strains 

A total of 70 Shigella strains isolated from 1997-2009 
were obtained from the Institute for Medical Research, 
Malaysia. However, only 40 strains S. sonnei from epide- 
miologically unrelated shigellosis cases were recovered 
from years 1997-2000 and 2007-2009 and analyzed. 
The strains were isolated from stools of patients in dif- 
ferent hospitals in Malaysia. Each strain was exclusively 
from one patient. The shigellosis cases were from differ- 
ent regions, although mainly from northern region of 
Peninsular Malaysia. 

Biotyping 

Biotyping was performed using fermentation of rham- 
nose and xylose and hydrolysis of ortho-nitrophenyl-/ 
I-D-galactopyranoside (ONPG), and biotypes were desig- 
nated according to methods described by Nastasi et al. 
[19]. Five biotypes are defined (biotypes a, d, e, f, and g) 
based on the combination of the positive/negative bio- 
chemical reactions. 



Antimicrobial susceptibility testing 

Antimicrobial susceptibility testing was performed using 
Kirby-Bauer disc diffusion method according to guide- 
lines of the Clinical and Laboratory Standards Institute 
[20]. E. coli ATCC 25922 was used as the quality-control 
strain. The antimicrobials commonly used for treatment 
of shigellosis were tested using antimicrobial-impregnated 
discs from Oxoid, Basingstoke, UK: ampicillin (A, 10 (ig), 
chloramphenicol (C, 30 (ig), streptomycin (S, 10 |ig), tetra- 
cycline (T, 30 (ig), ciprofloxacin (Cip, 5 (ig), kanamycin 
(K, 30 (ig), ceftriaxone (Cro, 30 (ig), nalidixic acid (N, 
30 (ig) and trimethoprim-sulfamethoxazole (Sx, 25 (ig). 



Pulsed-field gel electrophoresis (PFGE) 

All S. sonnei strains were analyzed with the restriction 
enzyme Xbal (Promega, Madison, WI, USA) using stan- 
dardized PulseNet PFGE protocol, with extended run 
time of 25 hours in a CHEF Mapper system (Bio-Rad, 
U.S.A.) [15]. XZ?^I-digested Salmonella serotype Braen- 
derup H9812 was used as the molecular weight standard. 
Computer-assisted analysis of the PFGE banding patterns 
was performed with BioNumerics software version 6.0 
(Applied Maths, Kortrijk, Belgium). Banding patterns were 
compared based on the Dice coefficient, and a dendro- 
gram based on the VFGE-Xbal profiles was constructed 
by unweighted pair group method with arithmetic 
averages (UPGMA) algorithm. A band position tolerance 
of 1.0% and optimization of 1.5% was applied during the 
comparison of banding patterns. 
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Multilocus variable-number tandem-repeat analysis 
(MLVA) 

Seven VNTR loci, SSI, SS3, SS6, SS9, SS10, SS11, and 
SSI 3 previously reported by Liang et al [18] were used. 
The primers and conditions are indicated in Table 1. 
Two multiplex polymerase chain reactions (mPCR) were 
carried out; mPCRl consisted of primers for SSI, SS3, 
SS6, and SS9 while mPCR2 consisted of primers for 
SS10, SS11, and SS13. Each 10 ul PCR mixture con- 
tained IX PCR buffer, 3 mM MgCl 2 , 0.2 mM of each 
dNTPs, 0.05 to 0.2 uM of each primer, 1 unit of GoTaq® 
Flexi DNA polymerase (Promega, Madison, WI, USA), 
and 1 ul of DNA template. Crude DNA template was 
prepared as follows: A pure culture of S. sonnei was plated 
on Luria-Bertani agar and incubated overnight at 37°C. A 
single colony was removed from the plate, suspended in 
100 ul of sterile deionized water, boiled for 5 min and im- 
mediately cooled on ice. After centrifugation at 12,100g 
for 3 min, the supernatant was transferred into a new tube 
for subsequent PCR analysis. PCR conditions were as fol- 
lows: initial denaturation at 94°C for 5 min; 30 cycles of 
denaturation at 94°C for 45 s, annealing at 55°C for 50 s 
and extension at 72°C for 60 s; and final extension at 72°C 
for 5 min. 

The resulting PCR products were diluted 20 times 
with sterile deionized water and separated on an ABI 
3130XL genetic analyzer using the size standard Gene 
Scan LIZ 500 (Applied Biosystems). Amplicon sizes were 
converted to copy numbers using BioNumerics software 
version 6.0 (Applied Maths, Kortrijk, Belgium). Each 
unique allelic string was designated a unique MLVA 
type. A dendrogram was constructed by UPGMA clus- 
tering based on categorical coefficient. Minimum span- 
ning tree algorithm was used to construct a minimum 
spanning tree (MST) to determine phylogenetic pattern. 
The allelic profiles of 173 strains from Taiwan, China, 
Indonesia, Vietnam, and Cambodia, previously reported 
by Chiou et al [16] and Liang et al [18] were extracted 



Table 1 MLVA primers, number of alleles and Nei's index 
of diversity (D) 



Locus* 


Dye Forward primer 
(5' to 30 


No. of 
alleles 


Nei's 
diversity 
index, D 


SSI 


VIC_TOCCAGTACACCTCACTCG 


13 


0.82 


SS3 


6-FAM_CTGGGAGATGAACAGGAGGA 


16 


0.92 


SS6 


N ED_GAGTCGCTAAACGOTGOT 


17 


0.90 


SS9 


P ET_CGC AATC AGCA A A ACA A AG A 


12 


0.88 


SS10 


6-FAM_ACGGTGGGCmCTCTACCT 


6 


0.75 


SS11 


vic_ctggtccgggaga™tcg 


8 


0.83 


SSI 3 


pet_agacgctggc™tgacgat 


3 


0.42 



* Liang et al. (2007). 



and compared to the present data in the MST analysis 
based on the seven tested loci. 

The polymorphism of each VNTR locus was repre- 
sented by Nei's diversity index, calculated as 1 - J] (al- 
lelic frequency) 2 . Simpsons index of diversity (D) and 
Non- Approximated Confidence Interval (C.I.N.A.) were 
calculated as described previously by Simpson [21] to 
compare the discriminatory powers of PFGE and MLVA. 

Results 

Biotypes of S. sonnei 

Two biotypes were identified among the S. sonnei 
strains, i.e. biotype a (ONPG +, rhamnose +, xylose -) 
and biotype g (ONPG +, rhamnose -, xylose -). Out of 
the 40 strains tested, 29 (73%) and 11 (27%) were of bio- 
types a and g, respectively. 

Antimicrobial susceptibility profiles 

All the 40 strains were sensitive to kanamycin, ceftriaxone 
and ciprofloxacin. The resistance rates of other antimicro- 
bials are as follows: streptomycin (67.5%), tetracycline, 
(40%), trimethoprim-sulfamethoxazole (37.5%), ampicillin 
(10%), chloramphenicol (10%) and nalidixic acid (2.5%). 
Nine strains (22.5%) were susceptible to all the antimi- 
crobials tested. Fifteen strains (37.5%) were multidrug- 
resistant (MDR, defined as resistance to three or more 
classes of antimicrobials) with resistotype SSxT being 
dominant (n=10), followed by ACST (n = 2), ACSx 
(n = 1), SSxTN (n = 1) and ACSTSx (n = 1). Resistotype 
SSxT was exclusively associated with 10 out of 11 S. sonnei 
biotype g strains. The remaining biotype g strain 
showed additional resistance towards nalidixic acid be- 
sides having a core resistance type of SSxT. On the con- 
trary, the biotype a strains have eight different resistotypes 
(Figure 1). 

Genotypes based on PFGE 

The genetic relatedness of the strains was determined by 
PFGE. XZ?^I-digested S. sonnei DNA generated 40 repro- 
ducible unique PFGE patterns, each with 15-24 DNA 
bands, with a Dice coefficient, F ranging from 0.70 to 
0.98. Based on 80% similarity, eight clusters (CI - CVIII) 
were observed (Figure 1). Most clusters contained 
strains from different geographical regions. Similarly, 
most clusters contained strains isolated in different years 
except cluster CVIII which contained two biotype a 
strains isolated in year 2000 with identical resistotype 
ACST. Cluster CI contained three biotype a strains from 
Kedah. Cluster CII contained two biotype a strains 
which were sensitive to all the antimicrobials tested. Ten 
of the 11 biotype g strains were dispersed in cluster Oil 
and cluster CIV. Cluster III contained five biotype g and 
two biotype a strains. One of the biotype a strain even 
shared 93.3% similarity with a biotype g strain. Another 
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Figure 1 Dendrogram of S. sonnei PFGE patterns generated by UPGMA clustering method using Dice coefficient. Key: all sensitive i 
ampicillin (A), chloramphenicol (C), streptomycin (S), tetracycline (J), and trimethoprim-sulfamethoxazole (Sx). 



five biotype g strains were grouped in cluster CIV with 
one biotype a strain. Cluster V contained three biotype a 
strains. Two strains in this cluster were closely related to 
each other (95.5% similarity) despite being isolated seven 
years apart from two different locations (Perak and 
Kedah). Cluster VI and cluster VII contained six and 
three biotype a strains, respectively. Seven biotype a 
and one biotype g strains were not assigned into any 
of these clusters. 

Genotyping by MLVA 

MLVA based on seven VNTR loci was performed to fur- 
ther characterize the S. sonnei strains. The number of 



alleles and diversity of each locus are listed in Table 1. 
Among the seven loci tested, SSI, SS3, SS6, SS9, and 
SS11 were more variable compared to SS10 and SS13. 
Locus SSI 3 showed the lowest diversity (D = 0.42) while 
SS3 was the most diverse (D = 0.92). From our result, 
the biotype a strains shared a single common allele at 
locus SSI 3 whereas two common alleles were detected 
at this locus for S. sonnei biotype g. The combination of 
the seven VNTR loci with high and low degree of vari- 
ability was sufficient for molecular subtyping of the 
strains. The 40 <S. sonnei strains were discriminated into 
40 different MLVA types. Two pairs of single-locus var- 
iants, SSM14 and SSM15, and SSM35 and SSM36 were 



Koh et al. BMC Infectious Diseases 2012, 12:122 
http://www.biomedcentral.eom/1 471-2334/1 2/1 22 



Page 5 of 8 



detected. The strains in each of the pair were isolated in 
the same year but from different regions. Two main 
clusters, MI and Mil were observed from the dendro- 
gram generated (Figure 2). Cluster MI consists of bio- 
type a strains while cluster Mil consists only biotype g 
strains. There was no clear demarcation of the strains 
isolated from different years and localities. From the 
MST, it is apparent that the S. sonnei in this study were 
heterogeneous as the tree branched out extensively. 
Most MLVA types have a distance of two to four loci 
from each other and the biotype a and g strains were dif- 
ferentiated into two main branches (Figure 3). The low 
variability locus SSI 3 established the 'trunk' of the MST, 
separating the biotype a and g strains while the six 
remaining highly variable loci further discriminated the 
strains in each branch, establishing the 'twig and leaves' 
of the tree. 

S. sonnei is one of the major etiologic agents of 
travel-associated diarrhoea [22]. Therefore our results 
were compared with MLVA profiles of 173 strains from 
five other Asian countries to determine if these strains 
were closely related. Analyses of MST (Figure 3) show 



that most of the Malaysian biotype g strains were 
related to the Taiwanese strains than those strains iso- 
lated from neighbouring countries (Indonesia, Vietnam, 
and Cambodia). Only one Malaysian biotype g strain 
was distantly related to the Vietnamese strains, sharing 
common alleles at two VNTR loci. The biotype a 
strains were also related to strains from Taiwan. In par- 
ticular, three Malaysian biotype a strains shared com- 
mon alleles at five VNTR loci with the Taiwanese 
strains. 

Comparison of discriminatory power of MLVA and PFGE 
in genotyping of S. sonnei 

PFGE and MLVA displayed excellent resolution in dis- 
criminating the Malaysian S. sonnei. Both methods were 
equally discriminative with a Simpsons index of diver- 
sity, D of 1.0 (C.I.N.A. 0.989-1.011). 

Discussion 

The predominant biotype identified in this study was 
biotype a (n = 29, 73%) and the remaining strains were 
of biotype g (n = 11, 27%). These two biotypes were the 
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Figure 2 Dendrogram generated from S. sonnei MLVA suptyping by the UPGMA clustering method using categorical coefficient. Key: 
all sensitive (#), ampicillin (A), chloramphenicol (C), streptomycin (S), tetracycline (T), and trimethoprim-sulfamethoxazole (Sx). 
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Figure 3 Minimum spanning tree for the 40 Malaysian Shigella sonnei strains and strains from five other countries. For Malaysian strains, 
each circle represents one strain with a unique MLVA profile. Each circle is noted with the MLVA type of the strain. MLVA types in yellow colour 
are Malaysian biotype g strains whereas those in green are Malaysian biotype a strains. For the profiles of strains from five other countries, the 
circle size is not proportional to the number of strains belonging to the MLVA profile but it is proportional to the number of outbreaks or 
unrelated episodes sharing the same profile. Circles in grey colour denote global clonal group A, orange denotes group B, and blue denotes 
group C. Additional information on these foreign strains is provided [see Additional file 1]. The distance between two MLVA types is numbered. A 
distance of one locus between two closest MLVA types is marked by a thick black line, a distance of two loci is marked by a blue line, and a 
distance of three loci is indicated by a thin red line. A distance of four loci is indicated by a broken line and a distance of five loci is indicated by 
a light broken line. 

V J 



most commonly reported biotypes in Australia from 
1990-2009 [23]. In many other countries, an increasing 
prevalence of S. sonnei biotype g was reported. Izumiya 
et al [14] reported the prevalence of S. sonnei biotype g 
among travel associated cases in Japan while almost all 
biotype a isolates were from patients with travel history 
to Southeast Asia. In Korea, S. sonnei isolated during 
1977-1986 were of biotype a, whereas isolates in 1991- 
2000 were of biotype g [24], Similarly, S. sonnei isolated 
from cases acquired in Ireland, Italy, United States, and 
a few African countries also suggested the increasing 



prevalence of S. sonnei biotype g since the 1990s 
[9,25,26]. More Malaysian strains however need to be 
analyzed to obtain a more accurate representation of the 
distribution and prevalence of S. sonnei biotypes in 
Malaysia. 

Overall, PFGE analysis in this study indicated that 
S. sonnei biotype a strains were genetically more diverse 
than biotype g strains. The non-rhamnose fermenting 
biotype g strains were given further attention as there is 
a spread of a pandemic clone of biotype g strains across 
different continents [9,10]. Ten out of 11 Malaysian 
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biotype g strains in this study exhibited similar charac- 
teristics of the pandemic clone that is biotype g, resist- 
ance pattern of SSxT, and similar but distinguishable 
PFGE patterns. A number of these biotype g strains were 
highly similar in PFGE banding patterns. The banding 
patterns, by visual comparison, were similar to some 
biotype g strains from Ireland and Italy [13]. These Ma- 
laysian strains could be phylogenetically linked to the 
pandemic biotype g clone. The biotype g strains were 
isolated from the period 1997 to 2000 and year 2008, in- 
dicating that these strains persisted in Malaysia for at 
least a decade. 

In this study, more than half of the strains were resist- 
ant towards streptomycin (67.5%) and all strains were 
susceptible to kanamycin. This is in agreement with an 
earlier study in Malaysia by Hoe et al [12] on Malaysian 
strains collected during the period 1997-2000. Resist- 
ance to trimethoprim- sulfamethoxazole and tetracycline 
are commonly reported in S. sonnei [5,14,24,27]. How- 
ever the resistance rates to trimethoprim-sulfamethoxa- 
zole and tetracycline in this study were relatively low at 
37.5% and 40%, respectively as compared to those 
observed in Taiwan, Thailand, Japan and Korea 
[5,14,24,27]. Most of the strains in this study remained 
susceptible to ceftriaxone, ciprofloxacin, ampicillin, 
chloramphenicol and nalidixic acid. These results were 
consistent with a study on S. sonnei from Northeast Ma- 
laysia [8]. The low resistance rate to nalidixic acid in 
Malaysia is in concordance with a report by Izumiya 
et al [14] where resistance to nalidixic acid was less fre- 
quent in S. sonnei originated from South-east Asia. 
MDR S. sonnei strains persisted throughout the study 
periods. Selective antibiotic pressure may lead to the 
persistence of MDR S. sonnei strains in Malaysia [12] as 
ampicillin and trimethoprim- sulfamethoxazole are used 
for the treatment of Shigella infection locally [28]. To 
the best of our knowledge, all the 40 infected patients 
recovered from their illness. 

PFGE and MLVA showed comparable discrimination 
in subtyping of S. sonnei. Both techniques further dif- 
ferentiated the S. sonnei biotype a and g strains. 
However MLVA was better at grouping the strains on 
the basis of bio types, and the overall percentage of 
similarity among strains was low when subtyping was 
done using MLVA. Furthermore, MLVA subtyping was 
more rapid and less laborious compared to PFGE. In- 
terpretation of result was less subjective and results 
were more readily comparable between laboratories. 
All these suggest that MLVA may be a suitable com- 
plement to PFGE or even an alternative for routine 
subtyping of S. sonnei. 

Most of the PFGE and MLVA clusters contained 
strains from multiple geographical locations, some clus- 
ters even contained strains isolated from distant parts of 



the country although they were epidemiologically unre- 
lated. This concurred with the report of Pichel et al. [3] 
where epidemiologically unrelated Argentinean S. sonnei 
from very distant geographical areas were clustered to- 
gether by PFGE. Epidemiologically unrelated MLVA sin- 
gle-locus variants, and strains isolated seven years apart 
at different locations yet with highly similar PFGE pat- 
tern were also observed. These observations where 
strains with no apparent epidemiological linkage were 
clustered together may be due to travel of individuals 
within Malaysia, long-persisted dissemination of differ- 
ent clones of S. sonnei throughout the country, person- 
to-person transmission of a particular strain over an 
extended period with minor genetic changes, or a com- 
bination of these events. 

Based on the MST analysis, S. sonnei in the present 
study were heterogeneous, and a number of strains in 
this study were related to S. sonnei from Taiwan and to a 
lesser extent, to those from neighbouring countries. 
Since S. sonnei is frequently found responsible for travel- 
associated diarrhoea, it is not uncommon that transmis- 
sion of S. sonnei across countries and even continents 
occurred easily. A study on the global distribution of 
S. sonnei clones divided the strains from 50 countries, 
including Malaysia into three major clonal groups [29]. 
The 173 strains from five countries that were used for 
comparison in this study belonged to clonal groups A 
(n = 148) and C (n = 17) which were globally spread and 
clonal group B (n = 8) which was found in Europe, Africa 
and Asia [29] (Figure 3). One hundred and fifty one of 
these strains originated from 10 shigellosis outbreaks 
(Ol - O10) in Taiwan between 1996 and 2004. Another 
22 strains were collected from nine epidemiologically 
unrelated episodes (El - E9) in Taiwan, China, Indo- 
nesia, Vietnam and Cambodia between 1998 to 2005 
[16,18]. Based onVNTR loci SSI and SS6, the Malaysian 
strains in the present study also belonged to these three 
major clonal groups. These observations suggest the 
S. sonnei strains circulating in Malaysia belonged to dif- 
ferent clones that were spread worldwide. Although our 
observation may not represent the precise clonal struc- 
ture of S. sonnei in Malaysia due to the limited sample 
size, the present study has nevertheless demonstrated 
the wide genetic diversity of S. sonnei circulating in this 
country. 

Conclusions 

Shigellosis is endemic in Malaysia. Our results suggested 
that S. sonnei circulating in Malaysia were heterogeneous 
and generally derived from different clones. MLVA based 
on seven selected VNTR loci was rapid, reproducible 
and highly discriminative and therefore may comple- 
ment PFGE for routine subtyping of S. sonnei. 
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Additional file 



Additional file 1: Information on the 173 5. sonnei strains from five 
other countries. Information provided are extracted from previous study 
by Chiou et al. (2006), Liang et al. (2007) and Filliol-Toutain et al. (201 1). 
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